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Abstract 

We investigated wliether development of resistance to a Bt crop in the presence of a natural enemy would be slower than 
without the natural enemy and whether biological control, in conjunction with a Bt crop, could effectively suppress the pest 
population. Additionally, we investigated whether insecticide-sprayed refuges of non-Bf crops would delay or accelerate 
resistance to the Bt crop. We used a system of Bt broccoli expressing Cryl Ac, a population of the pest Plutella xylostella with 
a low frequency of individuals resistant to Cryl Ac and the insecticide spinosad, and a natural enemy, Coleomegilla maculata, 
to conduct experiments over multiple generations. The results demonstrated that after 6 generations P. xylostella 
populations were very low in the treatment containing C. maculata and unsprayed non-Bt refuge plants. Furthermore, 
resistance to Bt plants evolved significantly slower in this treatment. In contrast, Bt plants with no refuge were completely 
defoliated in treatments without C maculata after 4-5 generations. In the treatment containing sprayed non-Bt refuge 
plants and C. maculata, the P. xylostella population was low, although the speed of resistance selection to Cryl Ac was 
significantly increased. These data demonstrate that natural enemies can delay resistance to Bt plants and have significant 
implications for integrated pest management (IPIVl) with Bt crops. 
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Introduction 

Tlie commercialization of plants expressing insecticidal crystal 
(Cry) proteins from Bacillus thuringiensis {Bt) for insect management 
has revolutionized agriculture and become a major tool for 
integrated pest management (IPM) programs [1-2]. In 2011, Bt 
crops were grown on nearly 70 million ha in 27 countries in 2012 
[3]. Bt crops have provided economic benefits to growers and 
reduced the use of other insecticides [1,4—6], suppressed pest 
populations on a regional basis [7-9], conserved natural enemies 
[10] and promoted biological control services in agricultural 
landscapes [6] . However, the development of insect resistance is a 
major threat to the sustainable use oi Bt crops [11-12]. 

Since Bt crops were first commercialized in 1996, there is 
evidence that three lepidopteran pests have evolved resistance to 
Bt crops in the open field [13-15] and one case of a coleopteran 
pest [16]. Resistance to Bt plants is a serious concern, but the 
relatively few number of cases is in stark contrast to many cases of 
resistance to conventional insecticides, which has occurred much 
more rapidly [17]. Commonly proposed reasons for the few 
confirmed cases of resistance to Bt plants are the high dose of Bt 
proteins expressed in plants and the use of refuges of non-Bt plants 
that can serve as a pool of Bt susceptible alleles in the population 
[18-19]. 



Another possible reason for the relatively few cases of resistance 
to Bt plants could be their safety to natural enemies that help 
suppress pest populations. Numerous studies have investigated the 
effects of Bt crops and Cry proteins on natural enemies (predators 
and parasitoids) in the laboratory and field [20-21]. A meta- 
analysis has confirmed the safety of Bt proteins, especially when 
compared to traditional insecticides [10]. When negative effects on 
natural enemies have been observed with Bt proteins, they appear 
to be due to the poor quality of the host and not the Cry protein 
[22], but see Desneux et al., 2010 [23]. The safety of several Bt 
proteins has been verified in tritrophic studies conducted with Bt- 
resistant or non-susceptible herbivores that avoided the problems 
of prey-quality in some previous studies [22]. Allowing iJi-resistant 
hosts to ingest Bt proteins and then feeding the hosts to natural 
enemies (both predators and parasitoids) has revealed no effects on 
the natural enemies [24-26]. However, some reports continue to 
suggest natural enemies may be harmed by Bt proteins [27], but 
these have been challenged [28]. 

The conservation of natural enemies by the use of Bt plants 
could also influence the development of resistance to i?/ crops. This 
question was first studied by Gould et al. [29] in their conceptual 
and mathematical models on tritrophic interactions of a plant, an 
herbivore and a natural enemy. Their simplest conclusion was 
that natural enemies that increase diflferential fitness between 
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susceptible and resistant phenotypes on host plants will accelerate 
resistance; those that decrease the differential will delay resistance. 
Johnson et al. [30-31] carried out controlled studies of a parasitoid 
and a pathogenic fungus that attack Heliothis virescens on Bt tobacco 
and concluded that the parasitoid would likely delay the 
development of resistance to transgenic tobacco, while the 
pathogen would likely promote the development of resistance. 
MallampaUi et al. [32] discovered that different prey species of a 
generalist predator had different effects on the development of 
resistance by Leptinotarsa decemlineata to Bt potato: the presence of 
on(; prey species delayed resistance while the other accelerated 
resistance. Heimpel et al. [33] reported that the form of egg 
mortality could influence the rate of resistance, but the importance 
of egg mortality depended on other ecological processes in the pest 
population. Other simulation models reported that natural 
enemies could slow insect resistance to Bt crops or Bt pesticides 
[34-36]. As this summary indicates, the literature contains 
suggestions that natural enemies could delay or accelerate 
resistance, depending on whether there is a differential impact 
on susceptible or resistant phenotypes. 

In the present study we used a unique system, composed of 
broccoli plants transformed to express Cry 1 Ac protein, a 
population oi Plutella xylostella, a global pest of crucifers [37] with 
a low frequency of resistant individuals to CrylAc and the 
insecticide spinosad, and the predaceous ladybird beede, Cokome- 
gilla maculata, to conduct a multigenerational study in the 
greenhouse. Our objectives were to determine: (1) if a natural 
enemy can delay the development of insect resistance to a Bt crop; 
and (2) if biological control in conjunction with Bt crops can 
effectively suppress the pest population. In addition, to simulate 
field-realistic conditions for both the predator and prey, we 
sprayed refuges with insecticide in some treatments, but not others, 
and observed those effects on the development of insect resistance. 

Results 

Population Density of P. xylostella 

The predator, presence of a refuge, and use of a spray on the 
refuge each influenced the population dynamics of P. xylostella per 
Bt plant over the 6 generations of the i^xpcrimc'nt. A repeated- 
m("asures ANOVA, with generation and treatm("nt as factors, 
yielded a significant effect for generations (i^^n,) = 77.101, 
P<0.001), for treatments (i^i = 31.788, P<0.001), and tiie 
interaction term generation*treatments (-F12.656 = 16.250, 
/'<0.001). During the l*" generation, few P. xylostella were found 
on Bt plants (Fig. lA), and there were no significant differences 
between treatments using one-way ANOVA (^4 = 0.2.58, 
P- 0.900). By the 2"*^ generation, there was an average of 7 P. 
xylostella larvae and pupae per Bt plant in the treatment with only 
Bt plants, but the number of P. xylostella was still about 1 per plant 
in the other treatments (^4 = 3.767, P= 0.026). Bt plants were 
completely defoliated and control failure was evident in the Bt 
plant-only treatment at the 3"^*^ generation when the number of 
P. xylostella had increased to 5 1 per Bt plant {F4 = 8.667, P= 0.001). 

At the 4* generation, the number of P. xylostella increased to 12 
per PZplant in treatment Bt+R (75% P^plants -1-25% non-P^ refuge 
plants) and 29 per Bt plant in Bt-HSR (73°/, Bt plants -h25% 
spinosad-sprayed non-Bt refuge plants), significantly higher than 
those in treatments Bt-HR-HCm (75% Bt plants -h25% non-Pi refuge 
plants -I- predator) and BtH-SRH-Cm (75% P/ plants 4-25% spinosad- 
sprayed non-Pi refuge plants -I- predator) (p4 = 21.294, P<0.001). 
There was still only about 1 P. xylostella per Bt plant in treatment 
Bt+R+Cm at the 5'*' and 6* generations, significantiy lower than 
other treatments (5* generation: p4 = 59.203, P<0.001; 6*^ 



generation: p4 = 61.164, P<0.001). Most importantiy, over the 6 
generations of the test, only treatment Bt-I-R+Cm maintained <2 
P. xylostella per Bt plant at each generation, suggesting the 
important role that the predator played in maintaining a low pest 
population. 

When spinosad was used (treatment Bt+SR+Cm), the pest 
population was also maintained at a low level, except for a flare-up 
in the 5* generation when the P. xylostella population was >40 per 
Bt plant. Without the use of the predator (Bt+R), the pest 
population gradually rose and peaked at 61 per plant by the 6* 
generation, compared to being maintained at about 1 for all 
generations when the predator was used (Bt+R+Cm). When the 
predator was replaced by the insecticide spinosad (Bt-I-SR), the pest 
population increased more rapidly and peaked at 107 per plant in 
the 6''' generation, again showing the strong and lasting benefit of 
the predator in maintaining a low pest population. 

Examining the pest population on the refuge plants in the cages 
provided another indication of the overall performance of the 
treatments. Using repeated measures analysis, there were signif- 
icant differences for generations (P3. 740 = 29.341, P<0.001), for 
treatments (^4 = 60.288, P<0.001), and the interaction term 
generation*treatments (Pi4.96o= 10.071, P<0.001). On refuge 
plants, the density of P. xylostella varied between 5-120 P. xylostella 
per plant between generations and treatments (Fig IB). Only the 
populations in Bt+SR+Cm, which combined Bt plants, spinosad 
and the predator, were consistently the lowest and did not exceed 
20 per plant in any generation. Use of the predator alone 
(Bt+R+Cm) maintained a relatively low and stable pest population 
over all 6 generations. In treatments of 100% non-Pi refuge plants, 
P. xylostella densities exceeded 100 per plant at each generation 
despite keeping only 3 defoliated plants with their larvae and 
pupae in each cage to reduce the overall population. Therefore, 
we did not include the treatment of only non-Pi plants in 
Figure IB. 

Population Density of C. maculata 

In treatments Bt+R+Cm and Bt+SR+Cm, predator populations 
generally increased as the pest, P. xylostella, populations increased. 
A repeated-measures ANOVA yielded a significant effect for 
generations (Pi.s^s = 1 1.873, P= 0.002), for treatments (Pi = 9.410, 
P= 0.022), and the interaction term generation*treatments 
(Pi 838 - 1 1-225, P= 0.003). Only a few C. maculata adults were 
found in the T' generation because only 3 pairs were released at 
the start of experiment (Fig. 2A). Predator populations on Pi plants 
in treatment Bt+R+Cm remained about 1 per plant because of the 
low pest population, especially on Pi plants through the 6* 
generation (Fig. lA). In treatment Bt+SR+Cm, predator popula- 
tions remained about 1 per plant until the 5* generation when 
they increased to >4.5 in the 5''' and 6* generations (Fig. 2A). 
Mean values for the 5th and 6th generations by the independent- 
test between the two treatments, respectively, differed significantiy 
(G5: i(6) = 4.562, P= 0.004; G6: i(6) = 6.268, P= 0.001). This likely 
maintained the pest population on the Pi plants at a low 
population (Fig. lA), although resistance increased (Table 1). 

Predator populations were generally higher on the refuge plants 
(Fig. 2B) than on the Pi plants (Fig. 2A), likely reflecting the higher 
pest density on these plants. For the three treatments Bt+R+Cm, 
Bt+SR+Cm, and non-Bt+Cm, there were no significant differenc- 
es in predator density in most generations. The repeated measures 
showed significant effects for generations (P3.460 = 12.234, 
P<0.001) and the interaction term generation*treatments 
(p6.92o = 3.241, P= 0.012), but not for the treatments (Pj^ 2.181, 
P= 0.175). 
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Figure 1 . Plutella xylostella populations on Bt plants (A) and on non-Bt refuge plants (B) in greenhouse cages (Means ± SEM). Bt+R: 
75% Bt and 25% non-fif refuge plants; Bt+R+Cm: 75% Bt and 25% non-Bf refuge plants with C. maculata; Bt+SR: 75% Bt and 25% spinosad-sprayed 
non-6f refuge plants; Bt+SR+Cm: 75% Bt and 25% spinosad-sprayed non-6f refuge plants with C maculata; Bt Only: 100% Bt plants only, Non-Bt+Cm: 
100% non-6f refuge plants with C. maculata. Different letters within the same generations denote significant differences (P<0.05, one-way ANOVA, 
Games-Howell test for G3 and G6 in A, G1, G3 and G5 in B, other comparsions by Tukey HSD test). 
doi:1 0.1 371/journal.pone.0090366.g001 



Resistance to CrylAc 

Survival of P. xylostella larvae indicates resistance and, on Bt 
broccoli, it was significantly affected when analyzed using 
repeated-measures for generations (F2.714 — 35.792, P<0.001), for 
treatments (i^6 = 81.199, P<0.001), and the interaction term 
generation*treatments (-fie.zsz ~ 6-719, P<0.001). Resistance lev- 
els reached 74.9% after only 3 generations in the treatment with Bt 
plants only, significantly greater than the survival in other 
treatments (^6 = 50.049, P<0.001) and >2x the next highest 
treatment of Bt-I-SR (75% 5i plants -I- 25% spinosad-sprayed non- 
Bt refuge plants) (Table 1). The high rate of resistance in the 



treatment with only Bt plants was sustained through the 6' 
generation. In treatment Bt-HSR, resistance also developed rapidly, 
with 56.()"/o P. xylostella surviving at the 4*'' generation, significandy 
greater than the treatments BtH-R, Bt-nRH-Cm, BtH-SRH-Cm and the 
control cages (Non-Bt Only and Non-Bt-I-Cm) (^6 = 33.319, 
P<0.001). The survival rates were 4-7% in treatment Bt4-R4-Cm 
(with predators) at all generations, but increased to 32.5% and 
39.5 in treatment Bt-HR (without j)rcdat()rs) at the 5* and 6* 
generations, respectively. There were significant differences 
among treatments (5* generation: i?6 = 26.302, P<0.001; 6* 
generation: ^6 = 38.850, P<0.001) with <10% survival in 
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Figure 2. Coleomegilla maculata population on Bt plants (A) and on non-Bt refuge plants (B) in greenliouse cages (Means ± SEM). 

Bt+R+Cm: 75% Bt and 25% non-6f refuge plants with C maculata; Bt+SR+Cm: 75% Bt and 25% spinosad-sprayed non-6f refuge plants with C 
maculata; Non-Bt+Cm: 100% non-6f refuge plants with C. maculata. Different letters within the same generations denote significant differences 
(P<0.05, one-way ANOVA, Tukey HSD test). 
doi:1 0.1 371/journal.pone.0090366.g002 



Bt+R+Cm, Non-Bt Only and Non-Bt+Cm, and >70% survival in 
Bt+SR, Bt+SR+Cm and Bt Only, while Bt+R liad nearly 40% 
survival in the 6* generation. These significant differences at the 
6* generation among Bt+R, Bt+R+Cm and Bt+SR+Cm highlight 
the ability of the predator to slow the rate of resistance to -S^ plants. 
Treatment Bt+R+Cm, which included the predator, had the 
lowest survival (6.0"/o) while the treatment without the predator 
(Bt+R) had 39.5% survival. Including the predator but using 
spinosad (Bt+SR+Cm), which has some toxicity to Cm [38], 
increased the development of resistance to Bt plants as indicated 
by the 72.8% survival (Table 1). 



Resistance to Spinosad 

In the treatments in which spinosad was used (Bt+SR and 
Bt+SR+Cm), only low levels of survival of P. xyhstella were 
detected. Survival of individuals removed from cages and fed 
spinosad-treated leaves was only about 1% at the 4* and 6''' 
generations. 

Discussion 

A large body of literature has evaluated the potential ecological 
risks of Bt crops to non-target organisms including natural enemies 
of insect pests [1,2,22,39]. However, as is the case with 
conventional insecticides, the interaction of natural enemies and 
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Table 1. Survival on Cry 1 Ac leaf (%) (IVleans ± SEM) of 
Plutella xylostella larvae from adults taken from cages. 



Treatments 


Generation 
















G3 


G4 




G5 




G6 




Bt+R 


3.3 ±1.86 A 


19.0 


±6.61 B 


32.5 


±14.55 BC 


39.5 


±14.06 B 


Bt+R+Cm 


4.7±3.18 A 


7.0±3.06 AB 


5.0± 


3.51 AB 


6.0 ±2.86 A 


Bt+SR 


29.3 ±4.01 B 


56.0 


±4.74 C 


87.5 


± 1 .44 D 


83.8 


±3.07 D 


Bt+SR+Cm 


3.3±1.93 A 


20.7 


±2.19 B 


36.8 


±6.85 C 


72.8 


±3.28 D 


Bt Only 


74.9 ±5.46 C 


72.0 


±6.76 C 


90.5 


±3.97 D 


93.5 


±1.19 D 


Non-Bt Only 


1.0±0.58 A 


1.3±0.33 A 


1.0± 


0.58 A 


3.7±1.33 A 


Non-Bt+Cm 


0.7±0.33 A 


1.0±0.58 A 


4.7± 


1.86 AB 


2.3 ±1.86 A 



Bt+R: 75% fif and 25% non-Sf refuge plants; Bt+R+Cm: 75% Bt and 25% non-Bt 
refuge plants with C. maculata; Bt+SR: 75% Bt and 25% spinosad-sprayed non- 
Bt refuge plants; Bt+SR+Cm: 75% Bt and 25% spinosad-sprayed non-6f refuge 
plants with C. maculata; Bt Only: only fif plants; Non-Bt Only: only non-fif plants; 
Non-Bt+Cm: only non-fif plants with C. maculata. 

Different letters within the same column denote significant differences (P<0.05, 

one-way ANOVA, Tukey's test). 

doi:1 0.1 371 /journal.pone.0090366.t001 

resistance development in their liosts has received far less 
attention. 

Our present results indicate that the predator, C. maculata, 
combined with non-Bt and unsprayed refuge plants, delayed 
resistance in the P. xylostella population to Bt broccoli plants, when 
compared to treatments without the predator. Some theoretical 
work predicted a similar conclusion, that natural enemies could 
slow the development of insect resistance to Bt crops or Bt sprayed 
on crops [30,34-36]. The present work, conducted over multiple 
generations, adds empirical evidence to the theory that natural 
enemies can slow the development of resistance to Bt crops. 

Arpaia et al. [35] reported that the distribution of C. maculata in 
potato fields is driven by the prey's density but, in spite of this 
aggregation, a significant number of them occur in areas where 
food density is low. In the present study, we observed some C. 
maculata adults and larvae on Bt plants, although the prey density 
was very low on diese plants. StiU the likelihood of a iJi-resistant P. 
xybstella individual reaching adulthood on transgenic plants will be 
lowered if there are natural enemies on Bt plants, and this, in turn, 
wiU slow the rate at which the resistance alleles increase in 
frequency. In our tests, C. maculata were found on Bt plants (Fig. 2B) 
where they likely encountered i?i-resistant P. xylostella and removed 
them, thus eliminating resistance alleles from the population. This 
is the likely cause of the slower rate of resistance in the treatment 
with the refuge and predator. Although spinosad is considered 
relatively non-toxic to natural enemies, it does have a toxic effect 
on C. maculata, although much less than "harder insecticides" [38]. 

Widespread planting oi Bt crops places strong selection pressure 
on pest populations, which could result in resistance to Bt and 
control failure [11,17-19]. Currently, an accepted cornerstone of 
resistance management programs for Bt crops in the United States, 
Australia, and elsewhere, is a refuge consisting of plants free oi Bt 
insecticidal proteins, thus allowing iJ/-susceptible alleles to persist 
in the population. When resistance is recessive, as has been the 
case so far for many key pests, then matings between resistant 
individuals emerging from Bt plants and those from refuge plants 
will greatly reduce the possibility of creating homozygous resistant 
offspring that can damage Bt plants [17]. In the United States, the 
Environmental Protection Agency (EPA) mandates noii-i?i' refuges 
of different sizes for different Bt crops. 



Empirical data have demonstrated the usefulness of refuges for 
Bt crops. Tang et al. [40] reported that pest population growth was 
influenced by refuge size, with the highest populations occurring in 
treatments that had either no refuge plants or aU refuge plants. 
That study also confirmed other models [41] that the development 
of resistance was inversely proportional to the size of the refuge. 
Most importantiy, Tang et al. [40] demonstrated that a balance 
could be struck by refuge sizes that would delay resistance while at 
the same time limit the growth of the pest population and damage 
to the crop, as has been observed in the field [7,9]. 

In the present study, results indicate that refuges slow resistance, 
but a sprayed refuge could accelerate it. These findings are in line 
with our previous field study [42]. In the treatment with only Bt 
plants without refuge plants, control failure and high insect 
densities were observed in the 3''' generation (two replications) and 
the 4* generation (the other two rephcations). Plants in the 
treatment of spinosad-sprayed non-i?i; broccoli plants (Bt-I-SR) were 
completely defohated between the 4* and 5* generations. 

In the treatment (Bt-I-SR-I-Gm) in which the predator and 
spinosad were used, 36.8% of the pest population survived on Bt 
plants by the 5* generation and by the 6''^ generation that had 
increased to 72.8% (Table 1). Spinosad killed most susceptible P. 
xylostella on refuge plants, and not enough susceptible individuals 
mated with the resistant individuals on Bt plants. Therefore, 
resistance development was quick although there were C. maculata 
in the treatment (Bt-I-SR-I-Cm). However, the pest populations 
were very low in this treatment with 16.5 P. xylostella larvae per 
refuge plant (Fig. lA) and 8.9 per Bt plant (Fig. IB) at the 6'*^ 
generation. In the 6**^ generation, we found > 1 00 C. maculata in 
each cage. This indicates that the high predator population 
controlled P. xylostella density to a low level, despite the fact that 
the pest had developed resistance to the Bt toxin. 

While farmers are concerned with reducing the likelihood of 
resistance, they are more immediately concerned with lowering 
the pest population to avoid crop injury. It is well worth noting 
that in our experiments we not only saw the lowest rate of 
resistance development in the prey when the predator was not 
decimated by the use of an insecticide, but also the lowest and least 
fluctuating pest population on Bt plants (Fig. lA) and low and 
stable pest populations on non-Bt plants in the refuge (Fig. IB). 
Thus, our data suggest that farmers can have sustainable 
management of pests if they combine Bt plants with biological 
control. 

In conclusion, this study provides empirical evidence to confirm 
the theory that natural enemies can delay resistance development 
to Bt plants, but also demonstrates that it can do so while 
maintaining a low pest density and low crop damage. Non-Bt 
refuges are necessary to delay resistance to Bt plants, but spraying 
refuges could accelerate resistance if sprays reduce the function of 
important biological control agents of the pest. We suggest that 
host-plant resistance with Bt plants and biological control can be 
fully compatible within an overall integrated pest management 
(IPM) program. Our results have significant implications for IPM 
and insect resistance management (IRM) for Bt crops. 

Methods 

Insects 

Three strains of Plutella xylostella were used to create a hybrid 
population for the cage tests: the susceptible Geneva 88 (SS), the 
Cry 1 Ac-resistant (CrylAc-RR), and the spinosad-resistant (Pearl- 
RR) [24,40,43] strains. The hybrid population was created by 
releasing 100 Fl RSI (G88 female x CiylAc-RR male), 100 RS2 
(G88 female x Pearl-RR male) and 300 G88 moths into a cage. 
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Table 2. Experimental treatments. 




Group 


Treatments 


Replications 


Bt+R 


75% Bt plants and 25% non-Bt refuge plants 


4 


Bt+R+Cm 


75% Bt plants and 25% non-Bt refuge plants and C. maculata 


4 


Bt+SR 


75% Bt plants and 25% non-Bt refuge plants treated with spinosad 


4 


Bt+SR+Cm 


75% Bt plants and 25% non-Bt refuge plants treated with spinosad and C. maculata 


4 


Bt Only 


100% Bt plants only 


4 


Non-Bt Only 


100% Non-Bt refuge plants 


3 


Non-Bt+Cm 


1 00% Non-Bt refuge plants and C. maculata 


3 


R: refuge. 

SR: refuge with spinosad. 
Cm: C. maculata. 

doi:10.1371/journal.pone.0090366.t002 



The total number of moths in the cage was 500 with a 1 : 1 ratio for 
female and male moths from each strain. Eggs were collected from 
the cage and put on artificial diet to rear Fl larvae. About 1,000 
moths from F1-F4 were used to produce a synthetic population. F5 
pupae were used in the selection experiments. The expected allele 
frequency of the synthetic population (square root of survival rate) 
was 0.1 for Cry 1 Ac and spinosad resistance. The mean survival of 
F5 larvae was 0.033 on CrylAc plants and 0.025 on spinosad- 
sprayed plants. Therefore, the actual initial allele frequency at the 
start of the experiment was estimated to be 0.057 (square root of 
0.033) for CrylAc resistance and 0.050 (square root of 0.025) for 
spinosad resistance [43]. While this is higher than would be 
expected in the field initially when Bt crops are released (perhaps 
10 ■ or lower [41]), we expected that these initial frequencies 
would allow us to see differences among the treatments in a 
reasonable time frame. 

Larvae and adult C. maculata were obtained from DuPont 
Pioneer (Johnston, lA) and maintained in a climatic chamber at 
27±1°C, 65±5% RH and 16:8 L: D at Cornell University's 
Department of Entomology at Geneva, NY. Both larvae and 
adults were reared on decapsulated eggs of brine shrimp, Artemia 
Jranciscana, (Brine Shrimp Direct, Ogden UT) and 1.5% agar 
solution provided separately as a water source. 

Transgenic broccoli plants and insecticide 

We used Brassica oleracea producing high levels of CrylAc for our 
Bt broccoli plants [44] . Bt broccoli plants with 8 true leaves were 
used, and analysis by ELISA indicated that the CrylAc protein 
level was 12. 33 ±1.62 Hg/g fresh leaf tissue (n= 7). To ensure the 
activity of the Bt broccoli, the plants were screened with P. xylostella 
neonates of Fl heterozygotes (G88 x CrylAc-RR) when plants 
were 4 to 5 wk old [43]. The CrylAc plants that killed 100% of 
neonates of Fl heterozygotes (SS x Cryl Ac-RR), indicating high 
levels of expression in the Bt plants, were used in the greenhouse 
experiments. Non-iJi broccoli (cv. Packman) was used in the 
refuge. 

A commercial formulation of spinosad (SpinTor 2 SC, 240 g 
[AI] /Liter) was used. Refuge plants in the treatments Bl and B2 
were sprayed in the cages using a small hand-held sprayer. The 
CrylAc plants were covered during the sprays to avoid drift. The 
concentration we used was 90 ppm, which was the lowest field 
dose hsted on the insecticide label. The insecticide-treated refuges 
were sprayed at the T', 3'^'' and 5* generations in order to keep the 
P. xylostella population in check in the cages. 



Experimental designs 

The selection experiment was conducted in greenhouses at 
Cornell University's New York State Agricultural Experiment 
Station and followed the general procedures used in previous 
studies [40,43]. Each cage was 1.8 m longxO.9 m widexl.7 m 
high and constructed of nylon netting. Adults, but not larvae, 
could easily move between the different broccoh types, which were 
separated by a nylon-netting barrier (0.9 m high). This arrange- 
ment simulated adjoining fields with frequent inter-field move- 
ment by adults but negligible movement of larvae. 

Seven treatments were included in the experiment (Table 2). 
There were 16 plants total in each cage (12 iJi plants plus 4 noa-Bt 
refuge plants). 250 F4 pupae of the synthetic P. xylostella population 
were released into each cage. Three pairs of 1-week old C. maculata 
were released into cages of the treatments Bt-I-R-I-Cm, Bt-l-SR-HCm 
and non-Bt-nCm when the P. xylostella larvae were 2""^ instars. 
During the pupal period, all old plants were cut at the base of the 
stem, and plants and pots (which might have pupae on them) were 
kept in the cages for a week in order to allow adults to emerge. 
New plants were introduced into the cage to provide foliage for 
egg laying or existing larvae that had defoliated plants. There were 
four replications (cages) of ffi-plant treatments (Bt-I-R, Bt-I-R-I-Cm, 
Bt-nSR, Bt-nSR-nCm and Bt Only) and three of the controls (Non- 
Bt Only and non-Bt-HCm). 

Data collection 

Older larvae (primarily 3rd or 4th instars) and pupae of P. 
xylostella, and pupae and adults of C. maculata on broccoli plants 
were counted every generation when larval and pupal densities 
peaked. To test for resistance, 30-40 larvae from non-Bt refuge 
plants were collected from each cage at the 3"*, 4* and 5* 
generations. The larvae were reared on diet in the laboratory, and 
then adults were allowed to mate and eggs were harvested and 
tested as described below. At the 6''^ generation, we collected at 
least 60 pupae on refuge plants from most of the cages. The 
survival of 2""* instars derived from the pupae collected in each 
cage was tested on CrylAc broccoh leaf disks in 30-ml plastic cups 
[43]. For each cage, a total of 100 larvae in 10 replications were 
tested on Bt broccoli, and non-Bt broccoli was used as a control. 
We also tested survival of 2"'' instars derived from Bl and B2 
treatments on spinosad-dipped broccoli leaf disks in 30-ml plastic 
cups at the diagnostic dose of 10 ppm [45]. A total of 100 larvae 
were tested (10 replications, 10 larvae/rep) for each cage. For both 
CrylAc plants and spinosad treatments, survival was determined 
after 3 days at 27± 1°C, 50± 10% RH and 16:8 L:D photoperiod. 
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Statistical analysis 

All statistical analyses were conducted using SPSS 17.0 
Windows [46]. Descriptive statistics are given as mean values 
and standard errors of the mean. Because the data fit the 
assumptions for parametric analysis, the repeated measures 
ANOVA was used with the factors of generations and treatments 
for analysis of the population of P. xylostella on Bt and non-Bt plants 
and the survival on Cry 1 Ac leaves of P. xylostella larvae from adults 
taken from cages. The survival rates were transformed by square 
root before analysis. For each generation, the differences of 
resistance development and the population of P. xylostella and C. 
maculata were analyzed by one-way ANOVA, and means were 
compared by Tukey HSD, if the data fit homoscedasticity, and 
Games-Howell test if not. Differences of the mean values in the 
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